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This paper presents the design of the modified sliding mode controller 
(MSMC) for the purpose of tracking the nonlinear system with mismatched 
disturbance. Provided that the performance of the designed controller 
depends on the value of control parameters, gravitational search algorithm 
(GSA), and particle swarm optimization (PSO) techniques are used to 
optimize these parameters in order to achieve a predefined system’s 
performance. In respect of system’s performance, it is evaluated based on the 
tracking error present between reference inputs transferred to the system and 
the system output. This is followed by verification of the efficiency of the 
designed controller in simulation environment under various values, with and 
without the inclusion of external disturbance. It can be seen from the 


simulation results that the MSMC with PSO exhibits a better performance in 
comparison to the performance of the similar controller with GSA in terms of 
output response and tracking error. 
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1. INTRODUCTION 

Electro-Hydraulic Actuator (EHA) system is one of the important drive systems in industrial sectors 
and most engineering practices. However, the system is highly nonlinear due to many factors, such as 
leakage, friction, and especially the expression of fluid flow through the servo valve [12]. With these 
attributes being present, the system’s performance will be reduced significantly. In order to overcome the 
issue previously highlighted, the controller utilized for the system should be robust enough to overcome the 
entire operating range. Therefore, a suitable controller needs to be designed for a positive performance to be 
shown by the electro-hydraulic actuator. 

The raised numbers of works conducted with electro-hydraulic actuator system have been proposed 
over the past decades ranged from linear control, nonlinear control to intelligent control strategies such as 
[3]-[6] have been proposed over the past decades. Based on the literature study in [7]-[10], sliding mode 
control (SMC) is found to be efficient and widely applied in nonlinear systems [11]. However, it has been 
noticed that most of the existing results regarding the sliding surface design are more emphasized on the 
matched uncertainties and disturbances attenuation, provided the insensitivity of the sliding motion of 
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traditional SMC to matched uncertainties and disturbance [12]. This indicates the presence of uncertainties 
and disturbances in the same channel as the channel of the control input. Despite this indication, the presence 
of these uncertainties in various practical systems may not be able to fulfil the well-known matching 
condition. In this study, the dynamic model and design requirement of electro-hydraulic actuator were 
obtained from National Institute for Aerospace Research, Romania [13]. As for the dynamic model used in 
this study, the track input disturbance took place on a different channel from the channel of the control input. 
It can be gathered that with mismatched system, not only the sliding motion of traditional SMC is severely 
affected by the mismatched disturbance, SMC is no longer capable of coping, in spite of its well-known 
robustness Given the significance of attenuating mismatched uncertainties and disturbances in practical 
applications, many authors have devoted themselves to designing the sliding surface for uncertain systems 
with mismatched disturbance. Additionally, a few number of literatures have used conventional SMC with 
some modification in its sliding surface[14]-[16]. These modifications are done to enhance the ability of the 
modified SMC (MSMC). 

It is also important to emphasise that estimation is seen as one way to improve the accuracy of 
MSMC. Therefore, the results which display the influence which will be posed by the parameters of the 
designed controller on the tracking performance will be shown in one of the sections in this paper. According 
to reviews, particle swarm optimization (PSO) is the most-applied technique for combination with sliding 
mode control in order to adjust its control parameters. Furthermore, gravitational search algorithm (GSA) is a 
relatively new nature-based optimization algorithm, which is in accordance to the law of gravity [17], [18]. 
However, both algorithm have not been specifically applied in the estimation of MSMC controller’s 
parameter for mismatched disturbance system, such as an electro-hydraulic actuator which has been used in 
this study. In fact, there has not been any researcher who has devoted themselves to conducting parameters 
estimation for MSMC for the improvement of its accuracy and performance. 

This paper presents an investigation into performance comparison between MSMC, which has been 
optimized using PSO and GSA. Sum Square Error (SSE) has been used as the objective function for 
obtaining the optimum value of the controller parameters. Essentially, a good tracking response will provide 
a small SSE value. Furthermore, comparative assessment of both optimization method for the system’s 
performance has been presented and discussed. Moreover, the main contents of this paper are summarized as 
follows: Section 2 illustrates the mathematical modelling of the developed system. Following that, the 
MSMC algorithm derivation is demonstrated in Section 3. Then, the discussion on optimization algorithms 
used in this study is presented in section 4. Observation results are discussed, compared, and presented in 
Section 5. Last but not least, summary and conclusions are presented in Section 7. 


2. RESEARCH METHOD 
2.1. EHA system modelling 

The actuator dynamic equation of electro-hydraulic actuator servo system is expressed as[13]. As a 
result, the differential equations which govern the dynamics of electro-hydraulic actuator servo system with 
external disturbance injected to its actuator is constructed as 


peak a) 


SS k f 3 FL 
X, = yg Mio ae ee + X3 ern 


(2) 


. Ss k c |Pa7x 
hy = 5X3 — Xa +7 [as k, (3) 
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The system equations can be expressed as a product of matrix such as: 


; 0 1 0 
4) [ik £ 0 
X2| = m m + ut |-— (4) 
° Cc Pa-x m 
X3 0 = Ss oe kı X3 Z, a as 0 
Ke Ke kc N 2 


where x, represents the displacement of the load while x, indicates load velocity and x3 represents the 
differential pressure p4 —pz2 between the cylinder chambers caused by the load. Meanwhile, F; represents the 
external disturbance inflicted on the system and it is possible for it to be a constant or a time varying signal. 
This system is a mismatched system due to the disturbance it inflicts on different channels from the control 
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input u .Table 1 display the parameters of electro hydraulic actuator servo system, which are represented 
by (1), (2) and (3). 


Table 1. Parameter of EHA Servo System 


Parameters Value Unit 
Load at the EHA rod (m) 0.33 Ns? /cm 
Piston Area (S) 10 cm? 
Coefficient of viscous friction (f) 27.5 Ns/cm 
Coefficient of aerodynamic elastic force (k) 1000 N/cm 
Valve port width (w) 0.05 cm 
Supply pressure (P,) 2100 N/cm? 
Coefficient of volumetric flow of the valve port (cg) 0.63 — 
Coefficient of internal leakage (k,) 2.38 x 10-3 cm>/Ns 
Coefficient of servo valve (k,) 0.017 cm/V 
Coefficient involving bulk modulus and EHA volume (k,) 2.5 x 107+ cm>/N 
Oil density (p) 8.87 x 1077 Ns? /cm* 


2.2. Controller design 

The study of sliding mode control has gained popularity in recent years as a methodology for 
controlling nonlinear systems with modelling uncertainties and external disturbances. It is known that the 
crucial and the most important step of SMC design is the construction of sliding surface S(t) which is 
expected to response desired control specifications and performance. The trajectories are enforced to lie on 
the sliding surface. Given the desired position trajectory as xq = [X1q,%2q.X3q]’ E R”, and defined 
žid =X2q, X2q =X3q. Vector of the system states are assumed measurable and defined as 


T ; y 
x= [x1,x2, x3] = [xp,%,P1 | € R”. The state error of the system is defined as 
ei = Xi — Xia (5) 


where i = 1 to 3 ande E R”. 

The control objective of this design a bounded control input u. Hence, the output position x, tracks 
as closely as possible the desired position trajectory xg in spite of various model uncertainties. In order to 
achieve the states of the system to track the desired trajectories at the same time, the function of a new sliding 
surface is proposed is similar as in [14]. 


SQ) =x, + x3 + Gei +c fe (6) 


where c,and c, are a strictly positive constant to be specified according to the desired dynamics of the 
closed-loop system, while e4 is dependent to c,and cz. The desired dynamic response for the system is given 
as S = $ = 0 when the sliding surface is moving. If S is forced to be zero, the desired dynamics is attained 
and the tracking error will converge to zero. Then, it can be obtained. 


6, = —X_ — X3 — C261 (7) 
To obtain the control law, the constant plus proportional reaching law method was applied as 


in [19] and [20]. The dynamics of the switching function are directly specified by this approach which is 
described by the reaching function of the form 


Š = —Qsign(S) — KS (8) 


where Q and K are a constant with positive value and sign(S) represents the signum function which have a 
piecewise function as below: 


1 ;S>0 
sign(S)=}0 ;S=0 (9) 
—1 ;S<0 
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Such that the derivative of (6) gives: 


SŠ = X, + £ + c16, + ce, (10) 
Substitute (5) in (10): 
$= Xz + X3 + c [tiia] + col xal (11) 


Substitute (1), (2) and (3) in (11): 


sS kı c |P,—X3 o- 
ke 2— k. Xz + a 5 k,u| + c,[%,-%1@] + 


21X41 —-X14] (12) 


The control law is obtained by substituting (8) in (12): 
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Then, the control law in (13) can be expressed as: 


U = Un + Ueq (14) 

where 
k f _s FL $ kı M `: : a 
u, = Ke | 2 prta” Woe we ee c1X1 +C1 Xia (15) 
EnA iamka C2X1 t C2X14 
and 
— Ke ae ae = 
Ueq = cky aa Qsgn(S) KS] (16) 


where up is nominal control for the nominal part of the system and Uegqis the equivalent control with 
switching function to ensure the robustness. 


2.3. Optimization algorithm 

In order to obtain a positive tracking performance, choice of control parameters for the designed 
controller is important. Therefore, the designed MSMC is integrated with PSO and GSA so that the optimum 
value of the control parameters can be attained. This is essential for any changes to be made on different 
values of external disturbance signal injected into the system. Figure | illustrates the block diagram of the 
proposed MSMC with PSO and GSA algorithm. 

MSMC, which is designed for electro-hydraulic actuator servo system, consists of four control 
parameters namely c,,c2,Q, and K. The value of these parameters should be chosen in order to reduce the 
amount of tracking errors. With this reason, for the improvement of the adaptation characteristic of the 
system, PSO and GSA are used to determine the optimum value of these parameters, specifically through the 
assimilation of modified sliding mode controller with these three algorithms. Besides, the way this task is 
performed by PSO and GSA is in accordance to SSE as an objective function. Moreover, a good tracking 
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response will result in a small SSE value. The formula of SSE is provided by: 
2 
SSE = Ek-1(Xk — Yrer ) (17) 
where k represents the number of iteration, while x, indicates the system output at k iteration and y,¢¢ is the 


reference input given to the system. Essentially, optimization algorithm is utilized with the purpose to 
minimize SSE value. 


Nonlinear 
system 


Figure 1. Block diagram of the MSMC with PSO, GSA and GPS algorithm 


2.3.1. Particle swarm optimization 

The particle swarm optimization (PSO) algorithm is a population-based search algorithm, which is 
in accordance to the simulation of the social behaviour of birds within a flock. It is a computational method 
which optimizes a problem through the effort of iteratively improving a candidate solution with regard to the 
given quality measure [21]. Individuals in a particle swarm emulate the success of neighbouring individuals 
and their own successes. Furthermore, this algorithm maintains a swarm of particles through the 
representation of potential solution done by each of the particles. The important parameters of PSO for this 
application are number of particle (i), number of dimensions which represent number of control 
variables (j), fitness, fit of global best position, (V(t)), fitness,the fit of personal best position (y;), fitness 
for each particle, (fit;(t)), particle’s position (x;), social coefficient (cz), cognitive coefficient (c1), inertia 
weight (w), lower bound value (Xmin), upper bound value (Xmax) and number of iteration (t). 


2.3.2. Gravitational search algorithm 

Gravitational Search Algorithm (GSA) is an optimization technique where its development is based 
on the concept of Newton theory and metaphor of gravitational kinematics. There are four specifications for 
each agent in this algorithm: position, inertial mass, and active and passive gravitational mass. The position 
of the mass corresponds to a solution of the problem. Its gravitational and inertial masses are determined by a 
fitness function. A particular solution is represented by each mass, and the algorithm is navigated through 
proper adjustment of the gravitational and inertial masses [17]. The position of the agent corresponds to a 
solution of the problem, and its mass is determined by using a fitness function Apart from that, the other 
main parameters for GSA are number of iteration (t), number of dimension which represents number of 
control parameters (d), and the best and worst fitness in population which is dependent on GSA’s objective 
function (represented by Best and Worst respectively) On the other hand, force (F4(t)), acceleration 
(a#(t)), velocity (v#(t)), position (x#(t)), mass (M,(t)), gravitational constant (G )), and Euclidean 
distance between two masses i and j, along (R;;(t)). 


3. RESULTS AND ANALYSIS 

Simulation experiment has been done using MATLAB/Simulink 2015 software. For the assessment 
of the controller’s response, a square-shaped wave signal is chosen as a reference trajectory for electro- 
hydraulic actuator system, as illustrated in [22]. The designed controller assists the system in tracking this 
signal so that smaller tracking error can be obtained. The external force, with a constant value of 
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F, = 10500N, is added as a perturbation to the system actuator. The implementation of optimization 
algorithms is conducted using parameters, number of particles, and the use of i, which represents the particles 
at 5, 10, 15, 20, and 25. The initial value of the number of iteration (t) is set at 10, followed by its increase to 
20, 30, 40, and 50 iterations. In case of the presentation of results in this section, the output plot is yielded by 
5 particles within 10 iterations, 15 particles within 30 iterations, and 25 particles within 50 iterations. These 
particles are chosen for PSO and GSA in order to observe the performance of the designed controllers, each 
with small, medium, and bigger number of particles and iterations. 

The conventional sliding mode control (CSMC) is generally applied to practical systems due to its 
simple design procedure. In this study, the designed CSMC was similar as the one used in [23]. Furthermore, 
it was integrated with PSO and GSA in order to allow the adjustment of the controller on these parameters to 
be conducted automatically. It was also with the purpose for the controller to be adapted to the dynamic 
changing of the system. The inclusion of CSMC in this study is to illustrate the difference between two types 
of sliding mode control in the same mismatched system, when result is concerned. As for the number of 
particles and iterations, these values were chosen in order to prove that the small number of these parameters 
will result in good tracking performance. The explanation regarding the acquirement of SSE from each 
combination of CSMC and MSMC with PSO and GSA respectively is displayed in Table 2. 


Table 2. SSE obtained from combination of MSMC with PSO and GSA 
CSMC MSMC 
PSO GSA PSO GSA 
i5,t10 232.1050 370.9705 32.4511 1609 
i15, t30 232.1031 308.0041 32.5014 746.6484 
i25, t50 232.1030 234.5374 32.7035 881.2971 


According to Figure 2, the combination of CSMC and MSMC with PSO and GSA enables the 
system output to track the reference input with different accuracy. At the beginning, the lowest number of 
particles and iteration are used which are 5 and 10. Initially, the lowest number of particles and iteration used 
are 5 and 10. As shown in Table 2, with this combination, MSMC-PSO exhibits the best system output in 
comparison to the system output of CSMC-PSO, CSMC-GSA, and MSMC-GSA, with SSE amounting to 
32.4511, 232.1050, 370.9705, and 1609 respectively. Even though MSMC-PSO produces oscillation at each 
corner of the output response, the values of SSE are the lowest. Besides, they are capable of tracking the 
reference input which is accurately provided with small error, in comparison to the capability of others. This 
indicates that PSO is able to comprehend the optimum combination of the parameter, even with the lowest 
SSE value. On the contrary, it is apparent that MSMC-GSA displays the most negative performance. 


— Ref 

—— CSMC-PSO 
——— CSMC-GSA 
—— MSMC-PSO | 7 
—— MSMC-GSA 


Displacement (cm) 


“0 l 2 3 4 5 6 7 8 9 10 
Time(sec) 


Figure 2. Position output for CSMC-PSO, CSMC-GSA, MSMC-PSO and MSMC-GSA with 5 particles 
within 10 iterations for external disturbance, 10500N 
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Figure 3 shows the position output for CSMC-PSO, CSMC-GSA, MSMC-PSO, and MSMC-GSA, 
with 15 particles within 30 iterations in the presence of external disturbance, which amounts to 10500N. 
MSMC-PSO has the tendency to track the provided reference input, with SSE value amounting to 32.5014. 
However, the oscillation of the output response still occurs at each corner of the position output responses. In 
comparison to the previous MSMC-PSO output response, no significant change occurs in this situation when 
the values of i and t are 5 and 10 respectively. Chattering occurs in the output response produced by CSMC- 
PSO and CSMC-GSA. This proves the existence of chattering when conventional sliding mode deals with 
any system with mismatched nonlinearities. Moreover, the amount of improvement shown by the SSE 
produced by CSMC-PSO and MSMC-PSO is not significant, despite the increase of the number of particle 
sand iterations from 5 and 10 to 15 and 30 respectively. 

The output response for CSMC and MSMC when the number of particles and iterations is increased 
to 25 and 50 respectively is illustrated in Figure 4. As seen in Table 2, the decrease of the SSE value for 
MSMC-GSA occurs, which is from1609, which is significant, to 746. This is followed by its increase to 881. 
Despite this increase, these values are still unable to match the performance outputs which is displayed by 
MSMC-PSO. Besides, the SSE value for MSMC-PSO still amounts to approximately 32, while the SSE 
values for CSMC-PSO and CSMC-GSA are maintained around 232 and 234 respectively. Provided the 
emphasis placed on the optimum combination which results to the production of the lowest SSE values for 
CSMC-PSO, CSMC-GSA, MSMC-PSO, and MSMC-GSA, the optimum combination of the number of 
particles (i) and iterations (t) is for the best performance output which results in the lowest SSE, as shown in 
Table 3. Based on Table 3, the combination of MSMC and PSO leads to smaller SSE, in comparison to its 
combination with GSA. The efficiency of MSMC can be seen through the production of SSE, when CSMC is 
used as a comparison tool in order to measure the accuracy of MSMC in terms of its management of the 
mismatched system. 


— Ref — Ref 
—— CSMC-PSO — CSMC-PSO 
—— CSMC-.GSA ——CSMC-GSA 


poe 
à Ar —MSMC-PSO l I WWD = —MSMC-PSO 
Www VA —MSMC-GSA [V 
j | 


| ——MSMC-G8A 
| 
| 
08 
—_—— 


G 06 Ẹ 06 Nw 

$ 04 2 04 a 7 
02 02 H= — os f 
ca l 2 3 4 5 6 7 8 9 10 u l 2 3 4 5 6 7 8 9 10 

Time(sec) Time(sec) 
Figure 3. Position output for CSMC-PSO, CSMC- Figure 4. Position output for CSMC-PSO, CSMC- 

GSA, MSMC-PSO and MSMC-GSA with 15 GSA, MSMC-PSO and MSMC-GSA with 25 

particles within 30 iterations for external particles within 50 iterations for external disturbance, 


disturbance, 10500N 10500N 


Table 3. Optimum combination number of particles,i and iterations, t which produced the lowest SSE 


CSMC MSMC 
PSO GSA PSO GSA 
Number of particles, i 25 25 25 10 
Number of iteration, t 10 50 20 30 
SSE 232.1030 234.5374 32.1762 743.2693 


4. CONCLUSION 

Based on the results and analysis discussed, the integration of MSMC and PSO produces better 
tracking performance with small tracking error, in comparison to the combination of the controller similar to 
GSA’s controller and the increasing number of iterations and agents used. With the presence of large external 
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disturbance, the integration of MSMC and PSO results in good performance. Besides, this integration results 
have high accuracy in tracking the reference signal provided to the system. Last but not least, for the next 
phase of this study, future work which investigates into the techniques required for improvement of the 
performance output of the proposed control strategy will be put into consideration. 
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